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ABSTRACT 

This work deals with the correlation functions and their 
methods of obtention. A method without a "pure" time delay 
is considered which uses Laguerre functions type filters. 

The possibilities of utilization of the nonsymmetric Laguerre 
type filter are analyzed. A very inexpensive set of circuits 
which are useful to obtain the mentioned filters for low 
frequencies signals are presented, together with their 
computer analysis and laboratory realizations. 
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I . INTRODUCTION 



Correlation analysis has become the most powerful tool 
to separate the signal from the noise, especially In those 
applications where the signals are hidden beneath a blanket 
of other signals and extraneous noise. 

The development of Instruments which can continuously 
sample even the most noisy signals and compute the correla- 
tion function simultaneously, have extended the use of 
correlation analysis Into radio astronomy, fluid and solid 
state physics, neurology, seismology, geophysical exploration 
and other areas [Ref. 3]* 

Radar and sonar, of primary importance in military appli- 
cations, use correlation detectors because It has been shown, 
for the case of the binary decision problem, that no matter 
what the philosophy employed for the design of a receiver 
is, the best detector is always a crosscomrelator detector 
[Ref. 5]. 

Theoretically, an Infinite signal to mcolse ratio can be 
obtained In the detection of periodic signals In noise 
through the correlation technique, but this would require 
an Infinite time of observation. In practice, by using 
crosscorrelation practical measurements can be made In 
reasonable time. 

Initially, theoretical considerations on correlation 
functions, their properties and methods to obtain them using 
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pure time delays will be presented. Next, this thesis will 
present a method which does not use a pure time delay and 
yields results comparable to classical methods. 

The mathematical formulation of the method and the cir- 
cuits that can be used are studied and compared with the 
standard methods. Finally an approximation to the necessary 
type of circuits was analyzed by digital computer methods 
and checked in the laboratory using actual components. 
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II. CORRELATION FUNCTIONS 



It is assumed that some physical process produces the 
time functions x^(t), X 2 (t), .... Xj^(t) simultaneously. 

It Is assumed further that the physical process is stationary, 
in hther words it is not changing over time. Also it is 
assumed that the time functions are not zero and do not have 
a d.c. component or they have zero mean if the time functions 
correspond to a stochastic process [Ref. 3]* 

A. AUTOCORRELATION FUNCTION 

The time average autocorrelation function for determinis- 
tic or known signals is defined for two different classes of 
signals, signals of finite average power (case II) and 
signals of finite energy or pulse signals (easel). 

For signals of case II the definition is: 

Oc(t) s lim ^ / X, (t) X, (t + t) dt (1) 

T-h» -t ^ 

and for signals of case I: 

(R,(t) = /”x (t) X (t + t) dt (2) 

>00 

If the time function under consideration is a stochastic 
process the autocorrelation function is defined as follows: 

R(t) = E[x^(t) x^(t + t)] (3) 
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If the stochastic process Is ergodlc 



= R(t) (il) 

According to the definitions given above. Figure 1 shows 
an elementary device to obtain the autocorrelation function. 
The process Implies multiplication of the signal by Its 
delayed replica and averaging. If the Instantaneous values 
of x^(t) and x^(t + t) are multiplied and tfae product, averaged 
over a sufficiently long time, the result is the autocorre- 
lation function If some means to vary t have been added to 
the circuit. 

1 . Properties of the Autocorrelation CRef. 3] 

The autocorrelation function has the following 
properties : 

a. Maximum. 

The autocorrelation function Is a maximum at 
X = 0, mathematically: 

1(R.(t)| <(il( 0) Schwartz Inequality (5) 

b . Value at x = 0 

The value of the autocorrelation at x = 0 is 

connected to the total power of the signal. If x^(t) is a 

2 

voltage, then for x = 0, x^^(t) is the average power of the 
signal (as measured on a 1 ohm resistor) . Hence(R(0) is 
the average power of signals of case II and(}^(0) is the 
energy in the pulse signals of case I. 
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c. Function 



The autocorrelation function Is a function of 
the form of which Is characteristic of the original 
gnal x^(t ) . 

d. Type of Function 

The value of the function for negative values of 
Ls Identical to that for the same positive values. This 
ans that the autocorrelation function Is an even function 
t: 



(R.(t) =(i^(-T) (6) 

e. Averaging Time » 

If the averaging time Is long compared with the 
clprocal of the lowest frequency In the original signal, 
en repeated measurements of the product for a given t will 
eld values very close to one another. 

f. Periodic Functions 

If x^(t) contains any periodic frequency compon- 
ts, then (ii( t) will contain each of these frequency compon- 
ts. So for certain time functlons,(R,(x) Is known If the 
Iglnal signal x^(t) Is given. 

g. Physical Processes 

Any physical realizable process that produces 
time function like x^(t) Is such that the value of the 
notion at time t becomes more Independent of the value at 
+ T as T gets larger. Thermal noise, and even the quantum 
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mechanical uncertainty principle, introduce randomness that 
cause a gradual loss of coherence between the signal and 
its delayed replica. This is true even for the oscillators 
in the most stable atomic clock. 

Hence for a signal arising from a real process 
the autocorrelation function approaches zero as t becomes 
sufficiently large . 

h. Coherence Time 

The value of t that causes a significant reduc- 
tion in the autocorrelation function is a measure of the 
coherence time of the original signal x^(t). 

i. The Frequency Domain 

A less obvious property of the autocorrelation 
function is indicated by: 

(R.(T)<^(t)(a)) (7) 

Equation (7) means that the autocorrelation 
function and the power spectral density are Fourier Transform 
pairs. Mathematically: 



4>(w) = /(K,(x)e dt 

— 00 


(8) 


(R.(t) ~ f <j)(a))e'^“^ do) 


(9) 



-.00 



The power spectral density for deterministic 
case II signals is: 
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( 10 ) 



4>(w) = 11m i X^^(o3) X^*(u)) 

- 11m ^ |X^^(o3)|^ 

T-*<o 

where the subindex T Indicates that the function Xj^(t) Is 
truncated and hence observed through a window of observation 
of duration T. 

The power spectral density for deterministic 
case I signals Is: 

(J)(o)) = X^(co) X^*(o)) = |X^(uj)|^ (11) 

The same property holds for stochastic processes, 
In which case; 



R.(t)^<()(u)) (12) 

and (J)(u)) = llm^E{|X (co)|^> (13) 

T-h» 

j . White Noise 

The autocorrelation function for very wideband, 

uniform (white) noise with r.m.s. value E Is an Impulse 

2 

function at T = 0 with amplitude E^. This means that one 
characteristic of such kinds of noise Is that Its Instantaneous 
value Is completely Independent of the value at any other 
Instant and that the coherence time of the process producing 
the noise Is practically zero or very short. 
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B. CROSSCORRELATION FUNCTION 



While the autocorrelation function of a signal is equi- 
valent to the traditional technique of power density spectrum 
analysis, there is no classical analogy for crosscorrelation 
analysis. Crosscorrelation is concerned with the relation- 
ship between two different signals that arise in some common 
process [Ref. 3]* 

The mathematical expression of crosscorrelation is: 

= 11m ^ / X, (t) x.(t + t) dt (14) 

T-«o -T 

1. Properties of the Crosscorrelation [Ref. 3] 

The properties of the crosscorrelation function are, 
in general, quite different from those of the autocorrelation 
function. 

a. Type of Function 

The crosscorrelation function is not an even 
function of T, which means: 

However, 

a relationship that has practical Importance in obtaining 
^"0^ negative delays. 

b. Coherence 

The crosscorrelation function approaches zero 
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as T approaches Infinity for signals that arise from real 
physical processes, due to the presence of noise and the 
uncertainty principle . 

Furthermore if x^(t) and X2(t) arise from two 
completely separate, unrelated processes 



c. Frequency Domain 

The same as in the case of autocorrelation 

functions : 




)^4>i2 



(w) 



(18) 



In this case, however, the physical meaning of 
clear. It can be called the spectrum of 
crosscorrelation of the time functions x^(t) and X2(t). 

Figure 2 shows some typical examples of cross- 
correlation between functions. 

Crosscorrelation analysis provides a powerful 
analytical tool. The ability to measure the degree to which 
the signals that arise from a common physical phenomenon 
resemble each other as a function of the delay time between 
them, can provide a much deeper Insight into the phenomenon 
being studied than a separate analysis of the properties 
of either signal alone. 
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III. APPLICATIONS OP CORRELATION 



A. DETECTION [REF. 4] 

Given two functions consisting of a signal S(t) and 
noise N(t ) : 



Xj_(t) = S^(t) + N^(t) 


(19) 


X2(t) = S2(t) + N2(t) 


(20) 



it is possible to crosscorrelate them with 1the following 
result : 



^ CS3_(t) + N^(t)][S2(t+T) + N2(t+T)] dt 









(t) +(2„ „ (T) +-® 



N1S2 



N1N2 



(T) 



( 21 ) 



All the terms of the above expression will be zero as 
't approaches infinity. 

Now it is assumed X^(t) = X 2 (t) and hence the four terms 
of the crosscorrelation become 



( 22 ) 

In the expression thus obtained (^ 22 ( t ) is the autocorre- 
lation function of the signal itself and is a function of 
T other than zero even for large values of t. The remaining 
terms approach zero as x becomes larger. So the first term 
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is a measure of the signal and the other ones are a measure 
of noise. 

If now it is assumed that the frequency of the given sig 
nal is known, an internal generator in the receiver can be 
used to generate the second input to the correlator and the 
two signals to crosscorrelate are: 



X^(t) = S^(t) + N(t) 


(23) 


X2(t) = S2(t) 


(2^0 


The above one is the typical case of a 
and if the two signals are crosscorrelated 


radar receiver, 
the result is: 




(25) 



Now the number of noise terms has been reduced to one 
and it is reduced to zero for large values of x. The pre- 
vious examples show how correlation techniques can be used 
to separate desired signals from the noise. 

B. ANALYSIS OP LINEAR SYSTEMS [REP. 3] 

It can be shown that it is possible to obtain a unit 
impulse response of a linear system by driving it with broad 
band-white-noise and crosscorrelating this input with the 
system output, as is shown in Pigure 3* The correlogram 
that results gives the same information about the system as 
if it were excited with an approximation to a delta function 
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and Its output recorded on an oscilloscope or x-y recorder. 

This technique has important practical applications 
because of its immunity to Internal system noise and because 
it permits to obtain the Impulse response of a linear system 
in the presence of a signal in the linear system. 

With real time methods of correlation it is possible 
using the method described above to keep the response of 
critical systems under constant surveillance and to make 
optimizing adjustments. 
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IV. METHODS TO OBTAIN CORRELATION FUNCTIONS 



There exist three methods to obtain correlation functions: 

A. METHOD THROUGH DEFINITION-TIME DOMAIN 

The easiest way to obtain a correlation seems to be the 
process indicated in its mathematical definition (1) and 
(14), this means: a delay Is given to one of the signals, 

then the other signal and the delayed one are multiplied and 
finally the product is averaged being the result^ ( t^) one 
point of the correlation function. Varying t it is possible 
to obtain a graph representing the correlation function. 

A basic block diagram to perform the process was shown 
in Figures 1 and 2. The essential feature of such devices 
is the temporal delay x between the two signals. When the 
signals are fluctuating voltages or currents, the provision 
of the "pure” delay is rather a problem because the available 
delay devices produce distortion on the signals for long 
delay times. 

Another problem of this method is that it is essentially 
discrete, and the number of points of the correlation function 
that can be obtained in a given time Interval depends upon 
the number of individual delay devices. 

Generally to avoid distortion, digital delay devices are 
used [Ref. 3] and hence it is possible to obtain a number 
of discrete values of the correlation function without 
distortion. 
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Although this method is not cheap it is widely used for 
receivers which use correlator detectors. 

B. INDIRECT METHOD-FREQUENCY DOMAIN 

Equations (8), (9), (10), (11) and (13) set up the 
relations between correlation functlons-tlme domain, and 
power spectral densitles-frequency domain. 

It is seen from the mentioned equations that correlation 
of two signals, or one signal with a shifted version of 
Itself, in the time domain is equivalent to a complex con- 
jugate multiplication of their linear spectra in the frequency 
domain. 

In other words, the autospectrum and the autocorrelation 
are related through the Fourier Transform, the correlation 
being the inverse transform of the power spectrum. The same 
relation applies to the cross-spectrum and crosscorrelation. 

Although this method is the faster one it is expensive 
because to obtain correlations , means have to be provided 
to obtain the power spectrums of the signals. Hence this 
method is typical of special signal processing machines 
which are designed to observe several characteristics of 
signals, and the power spectrums are not only used to obtain 
correlation functions. 

C. METHOD WITHOUT "PURE" DELAY-TIME DOMAIN 

This method has the advantage that a "pure" delay is 
not required. The correlation is obtained as a summation of 
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terms and the method is based upon the properties of 
orthogonal functions [Ref. 2], 

This is the only method essentially anaJLogous and as 
will be shovm in the following sections inexpensive 
correlators can be built for some special kind of signals. 
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V. CORRELATION WITHOUT A "PURE"' TIME DELAY 



A. MATHEMATICAL FORMULATION 
1. Analysis Problem 

It Is supposed that any correlation functlon^Rcx) , 
either autocorrelation or crosscorrelation, can be expanded 
in a series of orthogonal functions [Ref. 2], thus 

(R.(t) •= Z a 6 (t)[w(x)]^ 0 £ t < 00 (26) 

n=0 ^ ^ 

where the are polynomials which form an othonormal set 

with respect to the weight function w(x) in the range 
0 £ X < “• Then the coefficients, a^. In the series expansion 
can be found in the usual way and are given by 

= / (R,(t) e (x)[(o(x)]^ ^ dx (27) 

n 0 n 

Now the system of Figure 4 is considered. The 
linear network has an impulse response h^(4) . So the filter 
output is : 

00 

X (n;t) = / X, (t-u) h (u) du (28) 

. 0 

and the multiplier output is 

00 

x^(n;t) X 2 (t) = / x^(t-u) X 2 (t) h^(u) du (29) 
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being the output of the averaging device 
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x^(n;t) X 2 (t) = / x^TtMiy~xjrtT 5i^(u) du (30) 



0 



= / 0^ (u) h (u) du 

A n 




Now comparing (31) to (27) it is seen that the 



coefficient a^ is Identical with the average x^(n; t ) xJXD’ 
provided that the filter impulse response is so chosen that 



The filters whoise impulse responses satisfy Equation 



(32) are called "orthogonal filters," and if they could be 
built, the coefficients a^, would be obtained simply by 
switching in each of the filters in turn and determining the 
corresponding averages [Ref. 2]. 

Naturally, in any practical system only a finite 
number, say (N+1), of filters would be used.. 

2 . Synthesis Problem 

Equation (26) can be written 



h^(t) = 0^(t)[a)(t)]^"^ 



(32) 



00 

(K. (t) = [0)(t)]^‘^ Z a^ e^(T)[u)(T:)]^"^ (33) 



n=0 



or 




Z a h (t) 0 < t< 00 

-Ann — 



(3^0 
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In particular. If a finite number of terms are 
taken, say N, it can be written as an approximation to (33) 

<R„(t) = Z a h (t) (35) 

" n=0 " " . 

If Y is taken to be zero, (35) becomes 

®*«(t) = Co)(t)]“^ I a h (t) (35a) 

^ n=0 ^ ^ 

and if Y is taken to be 1/2 

= S a^ h^(t) (35b) 

n=0 

Figure 5 shows the complete synthesis of the method. 
It is possible to use half the number of linear networks 
calculating first the coefficients due to the particular 
signals to be correlated, and then after introducing the 
values of a^ the approximate correlation function can be 
seen on the screen of a CRT as a transient following the 
application of the input Impulse [Ref. 2]. 

A convenient system results if in fact the pulsing 
can be carried out repetitively at Intervals that are long 
compared with the filter time constants, and yet sufficiently 
frequently to give a steady trace on the CRT screen. This 
can always be arranged in practice [Ref. 2]. 

It can be seen from Equation (35a) that if Y is 
taken to be zero, the output of the circuit in Figure 5 has 
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to be multiplied by [o)(t)] ^ to obtain the correlation 
function. 

B. TWO REALIZABLE SYSTEMS 

The applicability of all the preceding theory depends 
on realizing a set of filters whose Impulse responses 
satisfy Equation (32). 

The polynomials and hence the Impulse responses 

h^(t) depend only on the weight function oj(t) and the 
parameter y. 

The weight function w(t) = a e~ was dhosen. 

1 . Case of Y = 0 

The transfer function H(s) of a sirngle R-C cell 
(low pass filter) is: 



■ T 



S + i 

T 



g 

s + a 



( 36 ) 



being 



a 




RC 



Hrc(s) corresponds in the time domain to the impulse 
response : 



hR^(t) = a e 



( 37 ) 



If now H(s) of a single C-R cell (high pass filter) 
is considered: 
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s 



( 38 ) 



s + a 



and the Impulse response In the time domain is: 



hcp(t) = -a e”“^ 



( 39 ) 



Now a number n of high pass filter cells are 
cascaded and between each cell a buffer amplifier is con- 
nected. Buffer amplifiers have ideally Infinite input 
impedance, zero output impedance and unity gain, and they 
Isolate the cells. 

Hence the transfer function in the s domain for the 
system of n high pass filter type cells is: 



to the input of the circuit considered above via a buffer 
amplifier, the resultant circuit is of the form shown in 
Figure 6 and its transfer function in the s domain is: 






(s + 



( 40 ) 



If a single low pass filter type cell is connected 




(s + a) 



( 41 ) 




( 42 ) 
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The inverse Laplace transform of the above one Is: 



a 






( 43 ) 



which can be written as: 



h (t) 
n 



a e 



-at 



n 

Z 

k=0 



nl 

(n-k) I 



C^at)^ 
k* kl 



(4H) 



= a e L (at) 
n 



where a e~ is the weight function and L^Cctt) are the 
Laguerre polynomials and the p'roduct is deTlned as Laguerre 
functions. 

The given expression for h^(^) satisfies Equation 
(32) and the conditions of orthogonality [Sef. 2] and hence 
the circuit of Figure 6 can be used to build the correlator. 

2. Case of Y = 1/2 

Another filter based on the same weight function is 
the lattice structure shown in Figure 7 • In this case the 
s domain transfer function is given by 



H^(s) 



a 

_2 

+ ^ 
2 



s - 



a 



n 






(US) 



and the impulse response is: 



h (t) = I 
n 2 




L (at) 
n 



( 46 ) 
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which is identical as the case of y = 0 except for a 
constant factor [Ref. 2]. 

3 . Frequency Domain Considerations 

The .nonsymmetric filter, case of Y = 0, was used with 
good results to build an autocorrelator [Ref. 2] which used 
N = 10. Hence N+1 = 11 filters were used, being the first 
one a low pass filter (n = 0) and the rest combinations of 
a low pass filter and a different number of high pass filters 
Isolated by cathode follower amplifiers. 

We saw in the ideal case, which ^ corresponds to 
perfect Isolation buffer amplifiers, that the transfer 
function in the frequency domain was (4l): 



H(s) 



__25_ (— s_) 
s + a s + a 



( 47 ) 



and because of the Isolation between cells it was taken as 
the product of the transfer function of a low pass filter 
and n R-C high pass filter cells. 

For the domain the above one becomes: 









, 2xn+l 
+ 03 ) 



( 48 ) 



being a function of a and n. 

Figures 8, 9, 10 and 11 show the transfer functions 
in the frequency domain of a R-C low pass filter, an R-C 
high pass filter, ten high pass R-C cells cascaded and 
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Ideally isolated, and of the ideal filters for n = 1 and 

■3 

n = 10, respectively. In all cases a = 10 was taken. 

The transfer functions were obtained with the aid of 
the CORNAP circuit analysis program [Ref. 1] for the case of 
individual cells and for the case of multiple cells cascaded 
taking the product of the individual transfer functions. 

From the curves it can be easily seen which is the 
frequency range for the incoming signals and which the band- 
width of the system. Also if the incoming signals are given 
it is possible to see which filters need postamplification. 

It has to be pointed out that the frequency domain 
considerations are only important in the determination of 
the coefficients (Figure 5) only case in which one of the 
input signals is convoluted with the Impulse responses of 
the filters. 

As the system is linear, by varying a it is possible 
to obtain filters for different ranges of frequencies. 

The correlator described in Ref. 2 used a = 0 type 
filters to autocorrelate a sine wave which frequency was 
152 Hz, being the basic time-constant of the filter a = 843 
sec~^. The results were practically Identical to those 
which can be obtained with conventional methods. 

C. A SET OF PRACTICAL FILTERS 

The buffer amplifiers needed to build the filters have 
to have infinite input resistance, zero output resistence 
and unity gain and the system mentioned in Ref. 2 used 
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cathode followers as buffer amplifiers. A way to avoid the 
use of such amplifiers was studied in .this work. 

1 . Factor of Isolation 

A filter of multiple cells can be constructed in 
such a way that the time constant is maintained the same 
through the cells and the resistance is multiplied by a 
factor of k from cell to cell as in Figure 12. It can be 
shown that this is equivalent to a certain amount of isola- 
tion between cells which depends on the value of the men- 
tioned factor: 




k = Factor of Isolation 



(49) 



The maximum value of k that can be used depends upon 
the necessary number of cells, and the values of the 
capacitors and resistors available. 

2. Computer Analysis 

In the Appendix a program is given which allows the 
calculation of the Laguerre functions for any values of a, 
n and t. Figures 13 through 23 show the Laguerre functions 

■5 

for a = 10 and for orders 0 up to 10. 

The CORNAP circuit analysis program was used to 
obtain the Impulse responses of circuits of several number 
of cells and different values of the factor of Isolation. 

Impulse responses of circuits up to seven cells (n=6) 
were obtained because the results had considerable errors 
for circuits with a large number of cells. 
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The results showed that circuits with a factor of 



Isolation of 10 have exactly the same Impulse response that 
the ideal ones. The lowest factor of isolation that gave 
reasonable approximation to the ideal impulse responses was 
foxmd to be 5. Hence, a factor of isolation of 10 can be 
used whenever it is possible and the use of values less than 
5 should be avoided. 

3. Laboratory Results 

Figures 24 and 25 show pictures of impulse responses 

■5 

of two kinds of circuits: first for a = 10"^ and a factor of 
Isolation of 10 for orders 0 up to 3» and second for the 
same a and k=5 for orders 0 up to 5 • 

It can be seen that the practical results obtained 
agree with the theoretical considerations and computer 
predictions. 

It was not possible to build circuits with a larger 
number of cells because of the characteristics of the avail- 
able pulse generator, oscilloscope and components. 
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VI . CONCLUSION 



The realizability of the method of Laguerre type filters 
or without a "pure" time delay to obtain correlation func- 
tions was shown together with a set of very simple filters 
to achieve them for low frequency signals. 

The following things require further investigation; 

Multipliers and adder that could be used with the filters 
presented in this thesis to build a correlator as the one 
shown in Figure 5* 

Factlblllty of building, through the techniques of Large 
Scale Integration, Laguerre type filters for several values 
of the base time constant alfa and with a behavior close to 
the ideal type of filter, which means that the circuits have 
to use buffer amplifiers. 
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APPENDIX 



PROGRAM TO COMPUTE THE VALUES 
OP LAGUERRE FUNCTIONS 



MAIN PROGRAM 

THIS PROGRAM CALCULATES THE LAGUERRE FUNCTIONS FOR 
ORDERS 0 UP TO N AS A FUNCTION CF THE ARGUMENT X. 

THE PRODUCT ALFA^TIME IS TAKEN AS THE ARGUMENT. 

BASIC TIME CONSTANT:K ALFA 

NUMBER OF DISCRETE VALUES OF THE ARGUMENT :M 

INCREMENT OF THE ARGUMENT : DELTA 

CIMENSICN Y(350) 

READ (5 t ICO) N 
READ(5, 100)M 
READ(5,101)DELTA 

STARTING VALUE OF ARGUMENT 

REA0(5,101)X 
READ(5, 102)KALFA 
L=N+1 

DO 3C I=1,M 
CALL LAP(Y,X,N) 

DAMP=KALFAv£XP (-XJ 
DO 20 J=1,L 
Y( J) =Y ( J J ' DAMP 
WRITE(6,2D0) Y( J) 

20 CONTINUE 

X=X+CELTA 
5C CONTINUE 

C 

100 F0RMAT(I2) 

101 FORMAT(F5.3) 

102 FORMAT! 17) 

200 F0RMAT(F9.2) 

C 

STOP 

END 
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SUBROUTINE 

SUBROUTINE LAP 
PURPOSE 

COMPUTE THE VALUES OF THE LAGUcRRE POLYNOMIALS 
L(N,X) FOR ARGUMENT VALUE X AND ORDERS 0 UP 70 
N 

USAGE 

CALL LAP(YtXiN) 

DESCRIPTION OF PARAMETERS 

Y -RESULT VECTOR OF DIMENSION N+1 CONTAINING 
THE VALUES OF LAGUERRE POLYNOMIALS OF 
ORDER 0 UP TO N FOR GIVEN ARGUMENT X 
VALUES ARE ORDERED FRCM LOW TC HIGH ORDER 
X -ARGUMENT CF LAGUERRE POLYNOMIAL 

N -ORDER OF LAGUERRE POLYNOMIAL 



C 



SUBROUTINE LAP(Y,X,N) 

DIMENSION Y(l) 

TEST OF ORDER 
Y(l)=l. 

IF(N) 1, 1,2 

1 RETURN 

2 Y(2)=l.-X 
IF(N-1)1,1,3 

INITIALIZATION 

3 T=l.+X 



DO 4 1=2, N 

4 Y(I + 11 = Y(I)-Y(I-1)+Y( I)-(T-^Y( I ) - Y ( I-l ) ) /F LCAT { I ) 
RETURN 
END 
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